Role of Cross Helicity in Cascade Processes of MHD turbulence 
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The cross helicity H =< v ■ b > characterizes the level 
of correlation between pulsations of the magnetic field b 
and the velocity field v. In the ideal three-dimensional 
magnetic hydrodynamics, it is an integral of motion along 
with the total energy E'^ = E'' + E\ E'" =< \v\'^/2 >, 
E'' =< |fep/2 >. The third integral of motion is the mag- 
netic helicity; but within the framework of this study, we 
consider the fields in which the average magnetic helicity 
is close to zero. 

As developed turbulence is random process one can 
expect that, if there are no special reasons, the devel- 
oped turbulence of conducting fluid (the MHD turbu- 
lence) should be characterized by a low level of cross- 
helicity. Exactly such a situation is usually considered. 
Interest in cross helicity arose after highly correlated pul- 
sations of velocity and magnetic field were found in the 
solar wind Analysis of the energy and helicity evo- 
lution in a freely decaying MHD turbulence showed that 
the helicity decays more slowly than the energy; hence, 
the degree of correlation of fields v and b determined by 



H/E' 



can increase m 



the correlation coefficient C 
time for the free decay [2]. 

By itself, MHD turbulence gives the possibility of de- 
veloping various scenarios. The specificity of the con- 
ducting fiuid hydrodynamics is the possibility of occur- 
rence of Alfven waves; it is assumed they play a key role 
in the turbulent cascade, which leads to the Iroshnikov- 
Kraichnan spectral law E''{k) ~ E^'ik) - fc-^/^ gQ. in 
the simulations on the grid 512'^ [5[, it was shown that 
no Kraichnan-Iroshnikov spectrum arises in the noncorre- 
lated turbulence without an external field, and the turbu- 
lence with the spectrum close to E''{k) ~ E'"{k) ~ 
is realized in the inertial interval. In Q the grid was 
expanded to 1024^ and the turbulence was considered 
both with and without the external field. For the MHD 
turbulence without the external field, the "—5/3" law 
was confirmed and a significant anisotropy for which the 
transverse pulsations follow the IroshnikovKraichnan law 
was revealed in the external field. 

The problem of the cross-helicity effect on the forced 
MHD turbulence was considered in [7| only in the context 
of the Alfven scenario (i.e., the turbulence that gives the 
"3/2" spectrum without the cross-helicity source). On 
the basis of the EDQNM approximation, it was shown 
that the system tends to the steady state in which the 
correlation coefficient proves to be much higher than the 
ratio of energy to cross-helicity input rates. In this case. 



the energy spectrum becomes steeper. 

The purpose of this work is to investigate the spectral 
properties of the developed isotropic (non- Alfven) MHD 
turbulence stationary excited by an external force, which 
also injects the cross helicity into the fiow simultaneously 
with the energy. 

We consider MHD turbulence with the magnetic 
Prandtl number of order of unity. On the integral scale 
L, an external force acts with a energy input rate equal 
to e. The same forces inject a certain cross-helicity input 
rate x ^-nd no magnetic helicity. We assume the equidis- 
tribution of kinetic and magnetic energy 5vf ~ 5bf « Ei 
within the limits of the inertial range. Also in the inertial 
range the energy flux is constant on any scale I and equal 
to the energy dissipation rate (the energy input rate) 



ti 



(1) 



where ti is the characteristic exchange time. In the the- 
ory of isotropic turbulence this time is usually estimated 
as the eddy turnover time r/ « I / 5vi . In the case of non- 
correlated velocity and magnetic-field pulsations, H — 
and this estimate can be accepted also for the MHD tur- 
bulence in which nonlinear interactions dominate instead 
of the Alfven waves. 

The basic idea of further arguments is that the injected 
cross helicity delays the spectral exchange (increases the 
time) 



ti 



dvi 



(2) 



and the delay coefficient ^/ is related to the correlation 
level for the velocity and magnetic-field pulsations on this 
scale. Hypothesis leads to the estimation of energy 
pulsations on the scale I in the form 



6vf 



2/3 



(3) 



In this case, the delay coefficient actually determines the 
deviation from the Kolmogorov "4/5" law 



(4) 



and Eq. ^ coincides with the Kolmogorov-Obukhov law 
Svf « (el)^/^ for a = 1. 

The delay of the cascade processes should lead to accu- 
mulation of the energy of turbulent oscilations (in com- 
parison with the noncorrelated turbulence for the same 
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power source). The application of estimate ([3]) to the 
energy-transfer scale I = L for which v\ Ki E gives 



2/3 



(5) 



The simplest assumption about the form of con- 
sists in the fact that, on the scales of action of external 
forces, the delay of the cascade transfer is determined 
by the quantity (1 — x/^); which is a characteristic of 
noncorrelation of perturbations introduced by external 
forces. Taking into account the quadraticity of the terms 
describing the processes of spectral transfer, we can as- 
sume that 'Cl ~ (1 ^ x/^)^^ which gives the estimate for 
the average energy of the stationary forced turbulence 
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Thus, the cross helicity blocks cascade energy transfer 
and leads to energy accumulation in the system. This ac- 
cumulation proceeds until the vortex intensification com- 
pensates the decreasing efficiency of nonlinear interac- 
tions. 

When assuming that the external forces inject the cross 
helicity in the turbulence with the given fiux x, it is neces- 
sary also to accept the hypothesis that the cross-helicity 
fiux is constant over the spectrum for a steady state, 
which gives 



Hi 



t 
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X, 



(7) 



at that, the cross-helicity-exchange time = ti^^ 
could not coincide with the energy-exchange time U. 
Let Ci be the correlation coefficient for the velocity and 
magnetic-field pulsations on the scale / 
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where the angular brackets mean averaging and Hi is 
the cross helicity on this scale. Thus, Hi « CiEi and the 
substitution in Eq. ([7]) gives 



CiSvf 
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(9) 



Comparing this expression with Eq. ([3|), we can relate 
the correlation coefficient to the characteristics of the ex- 
change rate on the corresponding scale 



(x) 



Gi 



(10) 



If the delay in the helicity and energy exchanges depends 
identically on the scale; i.e., Q ~ the velocity and 
magnetic-field correlation should be independent of the 
scale (and, on the contrary, the correlation independence 



of the scale means the identical dependence of coefficients 
and ^i on the scale). 
In the strongly correlated turbulence, the relation be- 
tween the exchange times and the eddy turnover time, 
which should be very large on the energy-transfer scale, 
decreases with the scale approaching unity on the dissi- 
pation scale. If it depends on the scale in the inertial 
interval by the power law 
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(11) 



the correction of the spectral energy distribution of pul- 
sations is unambiguously related to the parameter /i 



^2/3^2/3(1-1- A')^ 



(12) 



The conclusions about the role of cross-helicity in the 
stationary forced MHD turbulence are supported by re- 
sults of numerical simulations using a shell model of 
MHD turbulence Shell models describe the processes of 
spectral transfer in the developed turbulence with the 
help of a limited number of variables, each of which 
is a collective characteristic of pulsations amplitudes of 
the velocity field Un and the magnetic field i?„ in the 
wavenumber interval fc„ < |fc| < kn+i, where fc„ = A" 
and A is the interval (shell) width. The equations for 
collective variables are written to reproduce the "basic" 
properties of the initial equations of motion: the same 
kind of nonlinearity and integrals of motion. The shell 
models are an efficient tool for investigating the statis- 
tical properties of developed small-scale turbulence (see, 
for example, Q); in particular, they reproduce well the 
basic known properties of MHD turbulence and the small 
scale dynamo [9[. However, model [9] inherited the ba- 
sic disadvantage of cascade models associated with the 
method of describing the helicity (in these models, the 
different-sign helicity is attributed to shells with even or 
odd numbers n). In this work, we used a new model, 
which is obtained by a generalization on the MHD case 
of the model proposed in [l3| for the helical hydrody- 
namic turbulence. The model equations have the form 



dtUn = iK{K{U, U) - A„(B, B)) 



Re 



+ /„, (13) 



dtB„ = ^fc„(A„([/, B) - A„(B, t/)) 



Rm 



(14) 



where 



A„(A', y) — X'^{Xn+lYn+l + X*_^^Y*^i) - Xl'^_iYn 
+Xn^lYn-i + X*_iY*_i - X'^{Xl^^iYn + XnY^_^_i) 

the asterisk designates conjugation, while the super- 
scripts r, i are the real and imaginary parts. Without 
dissipation, the total energy E'^ = J2i\Un\'^ + |B„P)/2, 
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FIG. 1: Dependence of the average energy of the steadily 
excited MHD turbulence on the injected cross helicity x/e- 
In the inset, the energy evolution is shown at x = (the 
black line), x ~ 0.3 (the gray line), and x ~ 0.6 (the dotted 
line) . 

the cross helicity H = Y^iUnB^ + B„[/*)/2, and the 
magnetic helicity = EKHiK? " B^)/'^- ^ the 
magnetic field is zero, the hydrodynamic helicity = 
^k,iiiU*y — ?/,j)/2 is conserved. The distinctive fea- 
ture of the model is the possibility of the occurrence of 
arbitrary-sign helicity in an arbitrary wavenumber inter- 
val. In all simulations, the Reynolds number and the 
magnetic Reynolds number Re = Rm — 10^, and the 
shell width A = 1.618. Time is measured in dimension- 
less units equal to the eddy turnover time on the maximal 
scale. The force /„ operates only in the two highest shells 
(the larger scales) providing constant input of kinetic en- 
ergy e = 1 and cross helicity x- In Fig- 1; we show the 
time-average values of the total energy of the system ob- 
tained for various values of x/^ and agreeing well with 
estimate ^ to which the solid line in the figure corre- 
sponds. In the same figure, we show how the total energy 
E'^ of the system varies with time for various levels of 
cross helicity input rate. At x = 0, the time of attaining 
the quasi-steady state amounts to several vortex revolu- 
tions and, at X = 0.6, exceeds 100 dimensionless time 
units. In this case, both the average value of energy and 
the character of its oscillations vary. 

The energy accumulation is also accompanied by cross- 
helicity accumulation. In Fig. 2, we show the average 
values of integral correlation coefficient C =< H > / < 
E'^ >. It is substantial that the turbulence accumulates 
it at a low level of the injected cross helicity (x/s << 1) 
; i.e., the integral correlation coefficient greatly exceeds 
the ratio between the injected helicity and the injected 
energy. Thus, this tendency is inherent not only to the 
Alfven turbulence 7], but also to the isotropic (Kol- 
mogorov) MHD turbulence. At large values of x/s, the 
coefficient C tends to unity. Figure 3 shows how the en- 
ergy spectra vary with increasing the level of the cross 
helicity injected in the flow. We present the energy val- 
ues for each scale compensated on the quantity k^^^ . In 
such a representation, it is the horizontal line that corre- 



FIG. 2: Dependence of the average correlation level C —< 
H > / < E > on the injected cross-helicity level x/£- 
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FIG. 3: Compensated energy spectra. 

sponds to the spectrum k^^/^ . It can be seen that there 
is such a spectral range at x = 0, while both the en- 
ergy of each scale and the spectrum slope increase with 
X- It is of interest to trace directly the variation of the 
exchange time. In Fig. 4, we show the vortex turnover 
time and the energy and helicity-exchange times calcu- 
lated for each shell in the turbulence with a high level 
of cross helicity (x/e — 0.6). It is indicative that the 
exchange time on the integral scale exceeds the vortex 
turnover time by almost two orders of magnitude. This 
difference decreases with increasing the wave number and 
vanishes in the dissipative range. In the inertial range, 
the energy flux is also constant and the exchange time 
is unambiguously determined by the energy of pulsations 
of this scale; i.e., tn ~< m^j >. This means that the 
power law for i„ coincides with the slope for the energy 
of pulsations. In the case shown in Fig. 4, tn ~ ;0-89±o.02^ 
and the energy distribution in the inertial interval follows 
the law < ul ;J^-88±0 02 (the straight line in Fig. 3). 
The unexpected result is that the power law for the ex- 
change time is retained also in the dissipative interval 
(Fig. 4). The helicity-exchange time behaves similarly 
to the energy-exchange time, but it is always somewhat 
less. 

Thus, it is shown that the cross helicity blocks the 
spectral energy transfer in MHD turbulence and results 
in energy accumulation in the system. This accumulation 
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FIG. 4: Vortex-turnover and exchange times for the case 
of x/s = 0.6. The vertical Une corresponds to the inertial- 
interval 



proceeds until the vortex intensification compensates the 
decreasing efficiency of nonlinear interactions. The for- 
mula for estimating the average turbulence energy is ob- 
tained for the set ratio between the injected helicity and 
energy. It is remarkable that the turbulence accumulates 
the injected cross helicity at its low rate injection - the 
integral correlation coefficient significantly exceeds the 
ratio between the injected helicity and the energy. It is 
shown that the spectrum slope gradually increases from 



"5/3" to "2" with the cross heUcity level. 
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